There are binned and unbinned data that I would like to look at.

First the unbinned data. I want for each, number, biovol, flux, psd, psd with confidence intervals. read_csv(“dataOut/”)

library(tidyverse)
library(cowplot)
library(plotly)

Load all data

unbinned_DepthSummary <- read_csv("dataOut/unbinned_DepthSummary.csv")
Parsed with column specification:
cols(
  .default = col_double(),
  profile = col_character(),
  time = col_datetime(format = "")
)
See spec(...) for full column specifications.
unbinned_EachSize <- read_csv("dataOut/unbinned_EachSize.csv")
Parsed with column specification:
cols(
  profile = col_character(),
  time = col_datetime(format = ""),
  depth = col_double(),
  psd_gam = col_double(),
  vol = col_double(),
  sizeclass = col_character(),
  lb = col_double(),
  ub = col_double(),
  binsize = col_double(),
  TotalParticles = col_double(),
  nparticles = col_double(),
  n_nparticles = col_double(),
  biovolume = col_double(),
  speed = col_double(),
  flux = col_double(),
  flux_fit = col_double(),
  GamPredictTP = col_double()
)
binned_DepthSummary <- read_csv("dataOut/binned_DepthSummary.csv")
Parsed with column specification:
cols(
  .default = col_double(),
  profile = col_character(),
  time = col_datetime(format = "")
)
See spec(...) for full column specifications.
binned_EachSize <- read_csv("dataOut/binned_EachSize.csv")
Parsed with column specification:
cols(
  profile = col_character(),
  time = col_datetime(format = ""),
  depth = col_double(),
  psd_gam = col_double(),
  lb = col_double(),
  ub = col_double(),
  binsize = col_double(),
  vol = col_double(),
  TotalParticles = col_double(),
  nparticles = col_double(),
  n_nparticles = col_double(),
  biovolume = col_double(),
  speed = col_double(),
  flux = col_double(),
  flux_fit = col_double(),
  GamPredictTP = col_double()
)

Summary

Binned

plot_nparticles <- unbinned_DepthSummary %>% ggplot(aes(x = tot_nparticles, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1) + scale_y_reverse() + scale_x_log10() + guides(colour = FALSE) + labs(x = "#Particles/L") + theme_cowplot()

plot_biovolume <- unbinned_DepthSummary %>% ggplot(aes(x = tot_biovolume, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1) + scale_y_reverse() + scale_x_log10() + guides(colour = FALSE) + labs(x = "Unadjusted Biovolume (mass units/L)")+ theme_cowplot()

plot_flux <- unbinned_DepthSummary %>% ggplot(aes(x = tot_flux, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1) + scale_y_reverse() + scale_x_log10() + guides(colour = FALSE) + labs(x = "Unadjusted Flux (units/m^2/day)")+
   theme_cowplot()

plot_psd <- unbinned_DepthSummary %>% ggplot(aes(x = psd, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1, shape = 21, fill = "black", stroke = 1) + scale_y_reverse()  + guides(colour = FALSE) + labs(x = "Particle Size Distribution Slope") + theme_cowplot()

plot_speed <- unbinned_DepthSummary %>% ggplot(aes(x = tot_speed, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1) + scale_y_reverse() + scale_x_log10() + guides(colour = FALSE) + labs(x = "Unadjusted Flux Weighted Sinking Speed (m/d)") + theme_cowplot()


ggplotly(plot_nparticles)

ggplotly(plot_biovolume)

ggplotly(plot_flux)

ggplotly(plot_speed)

ggplotly(plot_psd)


cowplot::plot_grid(plot_nparticles,plot_biovolume, plot_flux, plot_speed, plot_psd) #%>% ggplotly()

plot_nparticles

plot_psd + scale_x_continuous(limits = c(-4.6, -2.5))

Binned

plot_nparticles <- binned_DepthSummary %>% ggplot(aes(x = tot_nparticles, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1) + scale_y_reverse(breaks = seq(from = 0, to = 2500, by = 100)) + scale_x_log10() + guides(colour = FALSE) + labs(x = "#Particles/L") + theme_cowplot()

plot_biovolume <- binned_DepthSummary %>% ggplot(aes(x = tot_biovolume, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1) + scale_y_reverse() + scale_x_log10() + guides(colour = FALSE) + labs(x = "Unadjusted Biovolume (mass units/L)")+ theme_cowplot()

plot_flux <- binned_DepthSummary %>% ggplot(aes(x = tot_flux, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1) + scale_y_reverse() + scale_x_log10() + guides(colour = FALSE) + labs(x = "Unadjusted Flux (units/m^2/day)")+
   theme_cowplot()

plot_psd <- binned_DepthSummary %>% ggplot(aes(x = psd, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1, shape = 21, fill = "black") + scale_y_reverse()  + guides(colour = FALSE) + labs(x = "Particle Size Distribution Slope") + theme_cowplot()

plot_speed <- binned_DepthSummary %>% ggplot(aes(x = tot_speed, y = depth, colour = profile)) + geom_point(alpha = 0.5, size = 1) + scale_y_reverse() + scale_x_log10() + guides(colour = FALSE) + labs(x = "Unadjusted Flux Weighted Sinking Speed (m/d)") + theme_cowplot()


ggplotly(plot_nparticles)

ggplotly(plot_biovolume)

ggplotly(plot_flux)

ggplotly(plot_speed)

ggplotly(plot_psd)


cowplot::plot_grid(plot_nparticles,plot_biovolume, plot_flux, plot_speed, plot_psd) #%>% ggplotly()

LS0tCnRpdGxlOiAiRXhwbG9yaW5nIHByb2Nlc3NlZCBVVlAtMjAxNyBkYXRhIgphdXRob3I6ICJKYWNvYiBDcmFtIgpkYXRlOiAiNC8yMC8yMDIwIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgpUaGVyZSBhcmUgYmlubmVkIGFuZCB1bmJpbm5lZCBkYXRhIHRoYXQgSSB3b3VsZCBsaWtlIHRvIGxvb2sgYXQuCgpGaXJzdCB0aGUgdW5iaW5uZWQgZGF0YS4gSSB3YW50IGZvciBlYWNoLCBudW1iZXIsIGJpb3ZvbCwgZmx1eCwgcHNkLApwc2Qgd2l0aCBjb25maWRlbmNlIGludGVydmFscy4KcmVhZF9jc3YoImRhdGFPdXQvIikKCmBgYHtyfQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShjb3dwbG90KQpsaWJyYXJ5KHBsb3RseSkKYGBgCgoKIyBMb2FkIGFsbCBkYXRhCmBgYHtyfQp1bmJpbm5lZF9EZXB0aFN1bW1hcnkgPC0gcmVhZF9jc3YoImRhdGFPdXQvdW5iaW5uZWRfRGVwdGhTdW1tYXJ5LmNzdiIpCnVuYmlubmVkX0VhY2hTaXplIDwtIHJlYWRfY3N2KCJkYXRhT3V0L3VuYmlubmVkX0VhY2hTaXplLmNzdiIpCmJpbm5lZF9EZXB0aFN1bW1hcnkgPC0gcmVhZF9jc3YoImRhdGFPdXQvYmlubmVkX0RlcHRoU3VtbWFyeS5jc3YiKQpiaW5uZWRfRWFjaFNpemUgPC0gcmVhZF9jc3YoImRhdGFPdXQvYmlubmVkX0VhY2hTaXplLmNzdiIpCmBgYAoKIyBTdW1tYXJ5CgojIyBCaW5uZWQKYGBge3J9CnBsb3RfbnBhcnRpY2xlcyA8LSB1bmJpbm5lZF9EZXB0aFN1bW1hcnkgJT4lIGdncGxvdChhZXMoeCA9IHRvdF9ucGFydGljbGVzLCB5ID0gZGVwdGgsIGNvbG91ciA9IHByb2ZpbGUpKSArIGdlb21fcG9pbnQoYWxwaGEgPSAwLjUsIHNpemUgPSAxKSArIHNjYWxlX3lfcmV2ZXJzZSgpICsgc2NhbGVfeF9sb2cxMCgpICsgZ3VpZGVzKGNvbG91ciA9IEZBTFNFKSArIGxhYnMoeCA9ICIjUGFydGljbGVzL0wiKSArIHRoZW1lX2Nvd3Bsb3QoKQoKcGxvdF9iaW92b2x1bWUgPC0gdW5iaW5uZWRfRGVwdGhTdW1tYXJ5ICU+JSBnZ3Bsb3QoYWVzKHggPSB0b3RfYmlvdm9sdW1lLCB5ID0gZGVwdGgsIGNvbG91ciA9IHByb2ZpbGUpKSArIGdlb21fcG9pbnQoYWxwaGEgPSAwLjUsIHNpemUgPSAxKSArIHNjYWxlX3lfcmV2ZXJzZSgpICsgc2NhbGVfeF9sb2cxMCgpICsgZ3VpZGVzKGNvbG91ciA9IEZBTFNFKSArIGxhYnMoeCA9ICJVbmFkanVzdGVkIEJpb3ZvbHVtZSAobWFzcyB1bml0cy9MKSIpKyB0aGVtZV9jb3dwbG90KCkKCnBsb3RfZmx1eCA8LSB1bmJpbm5lZF9EZXB0aFN1bW1hcnkgJT4lIGdncGxvdChhZXMoeCA9IHRvdF9mbHV4LCB5ID0gZGVwdGgsIGNvbG91ciA9IHByb2ZpbGUpKSArIGdlb21fcG9pbnQoYWxwaGEgPSAwLjUsIHNpemUgPSAxKSArIHNjYWxlX3lfcmV2ZXJzZSgpICsgc2NhbGVfeF9sb2cxMCgpICsgZ3VpZGVzKGNvbG91ciA9IEZBTFNFKSArIGxhYnMoeCA9ICJVbmFkanVzdGVkIEZsdXggKHVuaXRzL21eMi9kYXkpIikrCiAgIHRoZW1lX2Nvd3Bsb3QoKQoKcGxvdF9wc2QgPC0gdW5iaW5uZWRfRGVwdGhTdW1tYXJ5ICU+JSBnZ3Bsb3QoYWVzKHggPSBwc2QsIHkgPSBkZXB0aCwgY29sb3VyID0gcHJvZmlsZSkpICsgZ2VvbV9wb2ludChhbHBoYSA9IDAuNSwgc2l6ZSA9IDEsIHNoYXBlID0gMjEsIGZpbGwgPSAiYmxhY2siLCBzdHJva2UgPSAxKSArIHNjYWxlX3lfcmV2ZXJzZSgpICArIGd1aWRlcyhjb2xvdXIgPSBGQUxTRSkgKyBsYWJzKHggPSAiUGFydGljbGUgU2l6ZSBEaXN0cmlidXRpb24gU2xvcGUiKSArIHRoZW1lX2Nvd3Bsb3QoKQoKcGxvdF9zcGVlZCA8LSB1bmJpbm5lZF9EZXB0aFN1bW1hcnkgJT4lIGdncGxvdChhZXMoeCA9IHRvdF9zcGVlZCwgeSA9IGRlcHRoLCBjb2xvdXIgPSBwcm9maWxlKSkgKyBnZW9tX3BvaW50KGFscGhhID0gMC41LCBzaXplID0gMSkgKyBzY2FsZV95X3JldmVyc2UoKSArIHNjYWxlX3hfbG9nMTAoKSArIGd1aWRlcyhjb2xvdXIgPSBGQUxTRSkgKyBsYWJzKHggPSAiVW5hZGp1c3RlZCBGbHV4IFdlaWdodGVkIFNpbmtpbmcgU3BlZWQgKG0vZCkiKSArIHRoZW1lX2Nvd3Bsb3QoKQoKCmdncGxvdGx5KHBsb3RfbnBhcnRpY2xlcykKZ2dwbG90bHkocGxvdF9iaW92b2x1bWUpCmdncGxvdGx5KHBsb3RfZmx1eCkKZ2dwbG90bHkocGxvdF9zcGVlZCkKZ2dwbG90bHkocGxvdF9wc2QpCgpjb3dwbG90OjpwbG90X2dyaWQocGxvdF9ucGFydGljbGVzLHBsb3RfYmlvdm9sdW1lLCBwbG90X2ZsdXgsIHBsb3Rfc3BlZWQsIHBsb3RfcHNkKSAjJT4lIGdncGxvdGx5KCkKYGBgCgpgYGB7cn0KcGxvdF9ucGFydGljbGVzCnBsb3RfcHNkICsgc2NhbGVfeF9jb250aW51b3VzKGxpbWl0cyA9IGMoLTQuNiwgLTIuNSkpCmBgYAoKIyMgQmlubmVkCmBgYHtyfQpwbG90X25wYXJ0aWNsZXMgPC0gYmlubmVkX0RlcHRoU3VtbWFyeSAlPiUgZ2dwbG90KGFlcyh4ID0gdG90X25wYXJ0aWNsZXMsIHkgPSBkZXB0aCwgY29sb3VyID0gcHJvZmlsZSkpICsgZ2VvbV9wb2ludChhbHBoYSA9IDAuNSwgc2l6ZSA9IDEpICsgc2NhbGVfeV9yZXZlcnNlKGJyZWFrcyA9IHNlcShmcm9tID0gMCwgdG8gPSAyNTAwLCBieSA9IDEwMCkpICsgc2NhbGVfeF9sb2cxMCgpICsgZ3VpZGVzKGNvbG91ciA9IEZBTFNFKSArIGxhYnMoeCA9ICIjUGFydGljbGVzL0wiKSArIHRoZW1lX2Nvd3Bsb3QoKQoKcGxvdF9iaW92b2x1bWUgPC0gYmlubmVkX0RlcHRoU3VtbWFyeSAlPiUgZ2dwbG90KGFlcyh4ID0gdG90X2Jpb3ZvbHVtZSwgeSA9IGRlcHRoLCBjb2xvdXIgPSBwcm9maWxlKSkgKyBnZW9tX3BvaW50KGFscGhhID0gMC41LCBzaXplID0gMSkgKyBzY2FsZV95X3JldmVyc2UoKSArIHNjYWxlX3hfbG9nMTAoKSArIGd1aWRlcyhjb2xvdXIgPSBGQUxTRSkgKyBsYWJzKHggPSAiVW5hZGp1c3RlZCBCaW92b2x1bWUgKG1hc3MgdW5pdHMvTCkiKSsgdGhlbWVfY293cGxvdCgpCgpwbG90X2ZsdXggPC0gYmlubmVkX0RlcHRoU3VtbWFyeSAlPiUgZ2dwbG90KGFlcyh4ID0gdG90X2ZsdXgsIHkgPSBkZXB0aCwgY29sb3VyID0gcHJvZmlsZSkpICsgZ2VvbV9wb2ludChhbHBoYSA9IDAuNSwgc2l6ZSA9IDEpICsgc2NhbGVfeV9yZXZlcnNlKCkgKyBzY2FsZV94X2xvZzEwKCkgKyBndWlkZXMoY29sb3VyID0gRkFMU0UpICsgbGFicyh4ID0gIlVuYWRqdXN0ZWQgRmx1eCAodW5pdHMvbV4yL2RheSkiKSsKICAgdGhlbWVfY293cGxvdCgpCgpwbG90X3BzZCA8LSBiaW5uZWRfRGVwdGhTdW1tYXJ5ICU+JSBnZ3Bsb3QoYWVzKHggPSBwc2QsIHkgPSBkZXB0aCwgY29sb3VyID0gcHJvZmlsZSkpICsgZ2VvbV9wb2ludChhbHBoYSA9IDAuNSwgc2l6ZSA9IDEsIHNoYXBlID0gMjEsIGZpbGwgPSAiYmxhY2siKSArIHNjYWxlX3lfcmV2ZXJzZSgpICArIGd1aWRlcyhjb2xvdXIgPSBGQUxTRSkgKyBsYWJzKHggPSAiUGFydGljbGUgU2l6ZSBEaXN0cmlidXRpb24gU2xvcGUiKSArIHRoZW1lX2Nvd3Bsb3QoKQoKcGxvdF9zcGVlZCA8LSBiaW5uZWRfRGVwdGhTdW1tYXJ5ICU+JSBnZ3Bsb3QoYWVzKHggPSB0b3Rfc3BlZWQsIHkgPSBkZXB0aCwgY29sb3VyID0gcHJvZmlsZSkpICsgZ2VvbV9wb2ludChhbHBoYSA9IDAuNSwgc2l6ZSA9IDEpICsgc2NhbGVfeV9yZXZlcnNlKCkgKyBzY2FsZV94X2xvZzEwKCkgKyBndWlkZXMoY29sb3VyID0gRkFMU0UpICsgbGFicyh4ID0gIlVuYWRqdXN0ZWQgRmx1eCBXZWlnaHRlZCBTaW5raW5nIFNwZWVkIChtL2QpIikgKyB0aGVtZV9jb3dwbG90KCkKCgpnZ3Bsb3RseShwbG90X25wYXJ0aWNsZXMpCmdncGxvdGx5KHBsb3RfYmlvdm9sdW1lKQpnZ3Bsb3RseShwbG90X2ZsdXgpCmdncGxvdGx5KHBsb3Rfc3BlZWQpCmdncGxvdGx5KHBsb3RfcHNkKQoKY293cGxvdDo6cGxvdF9ncmlkKHBsb3RfbnBhcnRpY2xlcyxwbG90X2Jpb3ZvbHVtZSwgcGxvdF9mbHV4LCBwbG90X3NwZWVkLCBwbG90X3BzZCkgIyU+JSBnZ3Bsb3RseSgpCgpgYGAKCg==